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A multiband superconductor with a conventional phonon mechanism can develop an unconventional 
state with a nontrivial order parameter phase relation between the individual bands. We propose 
that such a state can explain recent experimental results on MgCNia , which are suggestive of both 
s-wave pairing symmetry and of unconventional superconductivity. We show that such a state gives 
rise to Andreev bound states, and to spontaneous currents, at surfaces and around impurities. We 
also investigate possible phase transitions between states with different order parameter symmetry. 
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\ Introduction. The discussion of unconventional su- 
perconductivity is usually limited to systems with p- or 

\ d-wave pairing. It has been noted recently that uncon- 
ventional properties can also appear in systems with a 
multicomponent "s-wave" order parameter (OP) arising 
from a conventional phonon mechanism Multi-band 

' superconductivity has been first studied shortly after the 
development of the BCS theory |^ and has experienced a 

' revival during recent years in connection with a variety of 
superconductors, including MgB2, RNi2B2C, Sr2Ru04, 
and many others 1). In a multiband superconductor, in 
which the intraband pairing interaction is most attractive 
in the s-wave channel, each Fermi surface (FS) sheet can 
develop a nearly-isotropic s-wave OP. If the interband 
Cooper pair scattering is repulsive, however, a nontriv- 
ial relative OP phase configuration can exist between the 

' bands §. In general such a state can break symmetries 
in addition to U{1), such as the time-reversal symmetry, 
responsible for unconventional behavior. 

We propose here that the antiperovskite compound ^ 
MgCNia could be a realization of this type of unconven- 
tional superconductor. The thermodynamic properties, 
and the observation of a Hebel-Slichter peak in NMR- 
measurements suggest conventional pairing. In con- 
trast, a zero-bias conductance peak has been observed in 
the quasiparticle tunneling spectra Interpreting this 
as evidence for subgap Andreev bound states at the sur- 
face would point towards unconventional superconduc- 

; tivity §,11 . These two seemingly contradictory observa- 
tions can be reconciled within a multiband s-wave picture 

\ as outlined above. Our aim is to consider this type of un- 
conventional superconductivity as a working hypothesis, 

\ and analyze observable properties that could be tested. 
Here we report briefly our main results; more details will 
be given later in a separate publication. 

In MgCNis, the important electron bands are mainly 
derived from a hybridization of the Ni Zd with the C 
2p orbitals. From band structure calculations the 
main features of the FS are three pill-shaped hole pock- 
ets centered around the X-points, and a nearly spherical 
electron FS around the F-point, see Fig. 0. Since these 
four FS sheets together comprise the largest part of the 




FIG. 1: Simplified view of the band structure of MgCNis. 

electron density of states (DOS) at the Fermi level, we 
ignore some minor additional sheets, namely the smaller 
cigar-shaped pockets along the T-R lines, and the "jungle 
gym" structure along M-R 

Nature of superconductivity in MgCNig . Throughout 
this letter we maintain the viewpoint that superconduc- 
tivity in MgCNis is mainly driven by the three symmetry- 
related bands around the X-points, since their Fermi ve- 
locity is quite low in some momentum directions, result- 
ing in a large DOS 1^. These bands furthermore consti- 
tute what distinguishes MgCNia from most "ordinary" 
BCS superconductors, and therefore are likely to contain 
the source of any unusual phenomena. Each of these 
three bands can individually develop an s-wave-like OP 
rja, identical in magnitude by symmetry, where a = 1,2,3 
is the band index. The repulsive Cooper pair scattering 
between these pockets, mediated by the Coulomb inter- 
action, stabilizes an unconventional state satisfying r]a = 
r]iuj~"^'^°'~^\ where m = ±1 and u) = e^'^*/^. The corre- 
sponding order parameter is fi„i = 7?i + cij™7]2 + ™'73. 

This can be understood if we view the Cooper pair 
scattering as an intrinsic Josephson coupling between dif- 
ferent bands, which gives rise to the Josephson energy 

E = J {C0S((/)12) -I- COS(023) + C0s((/)3i)} , (1) 

where (jjab is the phase difference between bands a and 
b. The coupling energies are identical by symmetry. If 
J is positive, this energy is minimized by 0i2 = (f>23 ~ 
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= ±27r/3. This two-fold degenerate state breaks time 
reversal symmetry. 

If the coupling between the X-pockets and the central 
F-band is strong enough, an unconventional state whose 
gap function has the fc-space structure 4'±i(k) ^ k^. + 
Lu^^ky + uj'^^kl can be induced on the F-band as well, 
even if this band on its own would favor standard s-wave 
pairing. The fj±i and '^±1 order parameters belong to 
the same two-dimensional representation Eg of the cubic 
point group Oh, but only "^±1 has all the features of a 
"usual" unconventional state, namely a non-trivial phase 
dependence of its gap function on the FS sheet, and even 
nodal points at k"^ = ky = k'^. On the X-bands, the 
(singular) nodal points would fall in the region between 
the FS sheets. 

Impurity bound states. The phase differences between 
the X-pockets, as well as the induced unconventional 
state on the F-band, can give rise to scattering-induced 
subgap states. Such states arise from an interference ef- 
fect of Bogoliubov quasiparticles that experience a phase 
shift when scattered between regions of the FSs with dif- 
ferent OP phase . For simplicity we will treat X-bands 
and F-band separately, neglecting quasiparticle scatter- 
ing processes between the two. In the presence of an 
impurity at the origin, the X-bands are described by the 
Bogoliubov-de Gennes (BdG) equations 



-E 

k'b 



9ah 






'9 ah 



(2) 

where a, 6 = 1 ... 3 label the three X-pockets, E^ de- 
notes the kinetic energy of electrons in band a with lat- 
tice momentum k, <I>J^ — (w^, u^) is a Nambu spinor, 
Qah — goSab + <7i(l ~ Sab) dcuotcs intrabaud and inter- 
band impurity scattering matrix elements, and N is the 
number of unit cells in the sample. In this approximation 
the gap functions — Ax(^°'~^ are spatially invariant, 
even in the vicinity of an impurity. Eq. (|^) is invariant 
under a 120° rotation about the [lll]-axis, since only 
the phase difference between the bands is physically rel- 
evant. This leads to the distinction of three sectors of 
eigenstates, corresponding to the three eigenvalues 1, oj, 
and w^, of the rotation operator, where only the first two 
contain subgap states. In each sector, Eq. (||) reduces to 
an effective one-band model, e.g. for the first eigenvalue: 
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where the symmetrized and "antisymmetrized" impurity 
scattering matrix elements are 



'^^9ab = 9q + 2gi, 5a = E ^9ab = 90 - 



91- 



(4) 



Note that the electron and hole components of the quasi- 
particles experience a different impurity scattering po- 



tential. For the second sector, gs and ga are inter- 
changed in Eq. (||). It is standard to solve these equa- 
tions for the bound state energies e, with the result 
e/Ax — ± cos {9s — Oa) , where the scattering phase 
shifts are given by 9s,a = a^'rcta.n [tt Nxiep) 9 s, a] , Nxiep) 
being the partial density of states of a single X-pocket 
at the Fermi level. On the F-sheet the analogous single- 
band formulation gives e/Ar = ±1/a/1 -|- (7rA^r(e_F)5r)^, 
where gr corresponds to the isotropic scattering strength. 
Hence strong s-wave impurity scattering can lead to sub- 
gap bound states, with energy close to zero, for all bands. 

Surface bound states Subgap states arise also from 
the reflection of quasiparticles at the surface. From the 
BdG equations, one obtains the bound state energies (in 
the limit of an infinite surface barrier) 

Im A* Ao,t 



±- 



(5) 



if the constraint min(|Ainp, |AoutP) > Re A*„Aout 
is satisfied If the order parameter ampli- 

tudes are equal, this expression reduces to e = 
±|Ain|cos[((/>in - 0out)/2] . Here Ain and Aout denote 
the order parameters "seen" by the incoming and out- 
going particle, respectively, and (pin and (^out are their 
respective phases. For specular scattering the momen- 
tum component parallel to the surface is conserved. At 
a [100]-surface, no subgap states appear on any band, 
since the phase difference is zero. At a [110]-surface, the 
scattering processes can connect two X-pockets with a 
phase difference of 27r/3, so that we expect subgap states 
with an energy e ~ Ax /2 here. On the F-band, with 
(i-like OP ^P+i given above, a straightforward calculation 
yields 



e - ±Ar [Sikjkpf ~ 1] /2, 



(6) 



where Ar is the maximum bulk gap, and fcz is the z- 
component of the quasiparticle momentum. For the [111] 
surface, the situation is qualitatively similar on the X- 
sheets, and for the F-band we get 



e = ±(Ar/2)sin2 6lcos3V', 



(7) 



where the angles 9 and ■0 provide a representation in 
spherical coordinates, so that 
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The above bound state spectra lead to a subgap DOS 



^ ' 3 I V 1 + 2A V 1-2A V2 y 
at a [110] surface, and 



7V(e) = ^arcosh( ^] of^-\ 



(9) 
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at a [111] surface. Here A = e/Ar, A = 
{A/ Ar){kpa/2Tr)'^ , and a is the lattice constant. Note 
that in addition to a sharp peak centered at e = 0, a 
broad range of subgap states with nonzero energy is gen- 
erated. 

Obviously imperfection on realistic surfaces would de- 
stroy the transverse translational symmetry and lead to 
diffuse scattering. Consequently, subgap states can ap- 
pear even at a [100] surface, and the bound state spec- 
trum would be generally more smeared. The zero-bias 
features observed in quasiparticle tunneling can be inter- 
preted as a signature of these Andreev bound states at 
the surface Q. In this way we have reconciled the ob- 
served standard form of the Hebel-Slichter peak, which 
results from the coherence effects of the s-wave pairing on 
each X-hand, with the presence of subgap surface states 
due to the unconventional phase structure of the OP. 

Additionally, we would like to emphasize that the pres- 
ence of subgap states is also connected with pair break- 
ing. In this sense, non-magnetic disorder would have an 
impact in this material on the the superconducting tran- 
sition temperature unlike in a conventional superconduc- 
tor. Interestingly, it was experimentally found that in- 
troducing C-vacancies lowers Tc, whereby this type of 
disorder is non-magnetic ] p^ . 

Ginzburg- Landau ( GL ) theory. A GL description al- 
lows a qualitative analysis of spontaneous currents, which 
are connected with the presence of subgap states [^ , and 
of possible phase transitions between states with differ- 
ent OP symmetry. In its general form, it requires three 
OPs 77a, a = 1 . . . 3, one for each X-pocket, and three 
OPs *o, *±i on the T sheet: 
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with m = 0, ±1, 



(11) 

where the same index m corresponds to phases of identi- 
cal symmetry on the F- and AT-sheets. We focus on the 
A-bands with the rj OPs, since their superconductivity 
is likely strongest. For repulsive Cooper pair scattering 
between different A-pockets, fjo has a lower transition 
temperature than the degenerate 77+1 and 77_i phases, 
and will therefore be neglected here. The GL free energy 
density is then 
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m=±l 



-I- 62|^+i|1?7-i| + K2 {(D+177+1)* • D_i^_i + c.c} , 

where D = {Dx,Dy,Dz) = V - iffA, and A is 
the vector potential. We also introduce the differen- 
tial operators V™ = {dx,(jj"^dy,uj~"^dz) and = 
{Dx,(jj"^Dy,uj~"^Dz). The current density is then 
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Currents around an impurity. If the impurity does 
not break any point group symmetries, it will introduce 
a term /imp = S {\fj+i [^ -I- |7?_i |^) S^{r) in the free energy. 
In an 77+1 bulk state also an ?7_i component will be in- 
duced around this defect. Denoting the deviations of 77+1 
and 77_i from their bulk value by dfi±, we arrive at the 
GL equations (to first order in S) 

[2b + ifiV^] Sr)^i + K^Vli Sfi+i = 0, (14) 
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These equations can be 



solved in Fourier space, giving rise to the currents 
4ie 5|77^fk|2^2 i^[Q*q] 



he {2a - Kiq^){2b - Kiq^) - K- 



(15) 



where q = {qx,ojqy,Uj'^qz) and Q = ql+ujqy + Lu'^ql. This 
corresponds to a angular structure 

j(r) ~ (x(y2 - z^),y[,^ „ z{x^ - y'))/r\ (16) 

which is an octupole. Consequently, the resulting mag- 
netic field has a basic quadrupolar structure with correc- 
tions of angular momentum I — 4, i.e. the field vanishes 
directly at the impurity site. 

Surface currents. At a boundary or interface, the fol- 
lowing surface terms enter the GL expansion: 

/surf = S{x') [Sl{\fl+f + + S2{rf+Ci,fl- + c.c.)] , 

(17) 

where x' = n ■ r, uj = n'^ + LoUy + Lo'^n'^, and n — 
(nxjUy^Uz) is the surface normal. A calculation similar 
to that for the impurity current yields a surface current 



j^|77^ri'lm[<:^*ii], 



(18) 
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which is localized on a length scale comparable to the 
correlation length. Here n = (rix,(jJny,uP'nz). A spon- 
taneous current is generated at a [110] surface, where 
Im[a)*n] ~ (—1,1,0), but not at [100] or [111] surfaces, 
where Im [cl;*n] is zero. 

Gompetition between d- and s-wave OPs. "^±1 is an 
externally induced OP that competes with the "natu- 
ral" ^'o OP on the F band, which has a different (s- 
wave like) symmetry than the inducing source. If the 
gap introduced by ^±1 is large enough, it will govern 
the quasiparticles on the F-Fermi surface entirely and be 
the only OP of this band. Weakening this induced pair- 
ing state, however, can open the door for a state with 
a finite ^'o-component. Therefore, at T = a quantum 
phase transition is possible for the F-band as a function of 
the impurity concentration. To illustrate this at least on 
a qualitative level, we consider the GL theory describing 
this competition, 

/r = E Kl*mP + 6:.|*™|^}+7l*+i|'l*oP 

m=0, + l 

+<?(*;i?7+i+c.c.), (19) 
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FIG. 2: Schematic phase diagram for the 'I'+i and "to order 
parameters on the central band, in the temperature vs. dis- 
order plane. Solid lines mark second-order phase transitions; 
the dashed line marks a region in which the transition may 
be first-order if the coupling between vp+i and ^'o is strong 
enough. 



where we regard 77+1 as an external OP that is assumed 
to be constant, and ignore any feedback of the F-band 
on the X-bands. For 77+1 = 0, only the s-wave state VPo 
is stable, and the coexistence of ^+1 and ^"0 is energeti- 
cally disfavored. Hence a'_^i and 7 are positive, while a'^ 
is negative at temperatures below Tco (< Tc), the transi- 
tion temperature of Vl/o- The phase of VPo is determined 
by additional terms involving which were omitted in 
the above expression. These are small, since VlZ-i is sup- 
pressed by the coupling to the X-bands, so that the phase 
of is only weakly pinned. Collective phase modes may 
hence be relevant near the transition jll]. 

We may influence the strength of the 5* +1 -component 
by non-magnetic impurities which with increasing con- 
centration suppress this unconventional pairing channel 
of the F-band, while the "s-wave" state is basically un- 
affected. The growing concentration of impurities, nimp, 
leads to the increase of a'j^^, causing the decrease of ^+1. 
A continuous onset of occurs when ag -I- 7|^'+ip = 
is satisfied. However, if the coupling coefficient 7, i.e. 
the repulsion between the two OPs, is sufficiently strong, 
then a first order transition is favored within this model, 
see Fig. |[ In both cases the new phase appearing cor- 
responds to the coexistence of the two OP, since ^+1 
enforced from "outside" can never be eliminated com- 
pletely. 

Conclusions. In summary, we have shown that the 
— at first glance contradictory — experimental results 
on MgCNis (Hebel-Slichter peak in NMR and zero-bias 
anomaly in tunneling spectra) can be reconciled by a 
multi-band unconventional superconducting state, where 
the OP on the dominant FS sheets is s-wave like, but 
with non-trivial phase relations between the bands. The 
resulting state violates time reversal symmetry. Subgap 



Andreev bound states and spontaneous currents are sig- 
natures of such a state that are accessible to experiment. 
Our specific predictions about the direction dependence 
of surface bound state energies can be tested in tun- 
neling experiments on single crystals. The spontaneous 
impurity-induced currents can be verified by zero-field 
muon relaxation measurements p^ , while the sponta- 
neous surface currents might be observed by scanning 
SQUID or Hall probes. Furthermore, the induced su- 
perconducting state on the central F-band can undergo 
a transition from a "d-wave" state to an "s -t- (i"-wave 
state with increasing impurity concentration. The first- 
or second-order transition (see Fig. ^) may or may not 
occur, depending on the precise values of microscopic pa- 
rameters. Such a transition may leave their signature in 
thermodynamic quantities such as specific heat, London 
penetration depth, or ultrasound absorption. Eventually 
we believe that MgCNia could be the exciting case of an 
unconventional superconducting state due to its multi- 
band structure 0]. 

We would like to thank D.F. Agterberg, Ch. Helm, 
H. Kusunose, and T.M. Rice for helpful discussions. This 
work was financially supported by the Swiss National Sci- 
ence Foundation. 



[1] 
[2] 

[3] 



[4] 
[5] 
[6] 
[7] 
[8] 



[9] 
[10] 

[11] 
[12] 



D.F. Agterberg et ai, Europhys. Lett. 48, 449 (1999) and 
Phys. Rev. B 60, 14868 (1999) 

H. Suhl et ai, Phys. Rev. Lett. 3, 552 (1959); J. Kondo, 
Prog. Theor. Phys. 29, 1 (1963); A. J. Leggett, Prog. 
Theor. Phys. 36, 901 (1966). 

H.J. Choi et ai, Nature 418, 758 (2002); M. Imada, J. 
Phys. Soc. Jpn. 70, 1218 (2001); S.-L. Drechsler et ai, 
Physica C 317-318, 117 (1999); D.F. Agterberg et ai, 
Phys. Rev. Lett. 78, 3374 (1997). 
T. He et al, Nature 411, 54 (2001). 
P.M. Singer et ai, Phys. Rev. Lett. 87, 257601 (2001). 
Z.Q. Mao et ai, preprint: cond- mat/0105280. 
H. Rosner et ai, Phys. Rev. Lett. 88, 027001 (2001). 
S.B. Dugdale and T. Jarlborg, Phys. Rev. B 64, 100508 
(2001); D.J. Singh and I.I. Mazin, Phys. Rev. B 64, 
140507 (2001); J.H. Shim, S.K. Kwon, and B.I. Min, 
Phys. Rev. B 64, 180510 (2001). 

S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 63, 1641 
(2000); T. Lofwander et ai, Supercond. Sci. Technol. 14, 
53 (2001). 

T.G. Amos et al, Sol. Stat. Comm. 121, 73 (2001); 

M.A. Hayward et ai, Sol. Stat. Comm. 119, 491 (2001); 

Z.A. Ren et ai, Physica C 371, 1 (2002). 

P. Kumar and P. Wolfle, Phys. Rev. Lett. 59, 1954 

(1987). 

G. Luke et al, Nature 394, 558 (1998). 



